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Abstract 

The Witten's supersymmetric interpretation of harmonic oscillator in one dimension 
(D — 1) is generalized to any D > 1. The emerging centrifugal barriers are regular- 
ized via a small imaginary shift of coordinates. The resulting VT symmetrized super- 
symmetry si (1/1) interrelates the spiked oscillator partners with different strengths 
of the barrier. The innovated form of the creation and annihilation operators is 
constructed. 
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1 Introduction 



Within the highly schematic formalism of the Witten's supersymmetric quantum 
mechanics |l|], many difficult methodical problems related to the existence and spon- 
taneous breakdown of the supersymmetric multiplets of bosons and fermions can be 
studied without any recourse to the less essential subtleties of the relativistic quantum 
field theory. The Witten's formalism recollects a years-long experience of mathemat- 
ical physicists with the ordinary differential Schrodinger equation in one dimension 
0. It enables us to make a choice between the harmonic oscillator V(x) = x 2 and 
other candidates. The former, "linear" model is exceptional as it renders possible 
the construction of the fermionic plus bosonic Fock space (cf. review M, p. 285). 

A significant drawback of the Witten's methodical laboratory has been described 
by Jevicki and Rodrigues |4j as a partial or even complete failure of the recipe when- 
ever the (super)potentials happen to be singular. We intend to re-consider critically 
the real depth of the latter difficulty, paying attention to the more dimensional har- 
monic oscillators which possess the "solvable" centrifugal singularity. 

Following the recent idea of Bender and Boettcher 0] and working within their 
so called VT symmetric quantum mechanics H [Tj, we shall replace our more- 
dimensional (and, admittedly, spiked) harmonic oscillator Hamiltonians by the spe- 
cific analytic continuation of their radial components ||. In the other words, we shall 
work with the complex, pseudo-Hermitian Hamiltonians [f| and assume for simplicity 
that their spectrum remains real and discrete. 

2 Complexification of coordinates 

There hardly exists a textbook on quantum physics which would not mention har- 
monic oscillator. The shape of its potential can be used as a fairly realistic model in 
atomic physics and quantum chemistry. Within nuclear theory, the exact solvability 
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of the A— body harmonic oscillator looks almost like a miracle and proves particularly 

helpful. In field theory, the simplicity of the annihilation and creation operators 

q + ip f q-ip 

a= — a 1 = — -=- (1) 

enables us to understand free fields as "local" excitations IT^a^ |0) of the vacuum. 

In the recent studies of quantum fields with broken parity [T(], |ll|] , various Hamil- 
tonians H have been considered and modified by the complex shift of coordinate, 

q — > r = r(x) = x — ie, x G (—00, 00). (2) 

The original symmetry with respect to the parity V was broken and replaced by a 
significantly weaker VT symmetry [H, VT] = where T mimicked the action of the 
time-reversal of fields. The Hermiticity of H was lost. 

Notably, once we work with a "quasi-parity" quantum number (3 = ±1/2, the 
usual spectrum of the Hermitian harmonic oscillator Hamiltonian remains the same 
for all its shifted, non- Hermitian descendants ||, 

E = } = An + 2(3 + 2, n = 0, 1, . . . (3) 

The preservation of the discrete, real and semi-bounded form of these energies can 
be attributed to the VT symmetry of the underlying Hamiltonian . Many exactly 
solvable VT symmetric models seem to support such a hypothesis JTT|, [12|] . 

Once we fix the shift e > in eq. (H), the one-dimensional harmonic-oscillator 
Schrodinger equation can be replaced by a family of its singular generalizations 



d 2 a 2 - I I A 
- — C (r) + (r) + r 2 C (r) = EC (r) (4) 

GfC it 

with the same boundary conditions. The parameter 

a = (D-2)/2 + £ + lu>0 

varies with the spatial dimension D and with the non-negative integer angular mo- 
mentum quantum number I = 0, 1, . . .. 
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Our previous energy formula (|3]) remains valid and continuous at a = 1/2 pro- 
vided only that we put (3 = ±a [|T3|]. Off the point a = 1/2 the spectrum of our new 
equation H^ip = Etp is richer than in the Hermitian case. The analytic, confluent 
hypergeometric origin of equation (Q) is to be seen in a new perspective. The ordi- 
nary Sturm-Liouville theory must be adapted to the new situation [I3|. The norms 
have to be replaced by the pseudo-norms [|nj. At the same time, our differential 
Schrodinger equation (|4]) stays safely regularized at any a > 0. 

In the generic case with a / 1/2, our VT regularization (e ^ 0) is only removable 
via the limiting transition e — > accompanied by a halving of the axis of coordinates. 
This means that we have to replace r(x) = x — ie by the radial and real r e (0, oo). 
We must cross out all the states with (3 < 0. The e — > return to the Hermitian 
cases D = 1 or D = 3 remains transparent. A few details are mentioned in Table 1. 
The weird difference between the linear and spherical Hermitian oscillators can be 
attributed to the mere physically well motivated difference in boundary conditions 
in the origin. 

This initiated our study. Its decisive encouragement came from the discovery 
of confluence of the different energy levels at certain couplings a. This anomaly re- 
mained unnoticed in ref. and proves closely connected to the (unavoided) crossing 
of the energy levels at certain strengths of the barrier |Tj|. One can easily verify that 
the latter puzzling phenomenon takes place at all the integers a, i.e., at all the even 
dimensions D — 2, 4, ... for the vanishing residual u = 0. 

Such an observation is slightly discomforting. In the language of linear algebra, 
the exceptional character of the integers a can be intuitively explained via the oc- 
currence of Jordan blocks in the canonical form of the non-Hermitian Hamiltonian 
matrix Jl6[| . In order to simplify the discussion, we shall skip this point completely 
and, everywhere in what follows, we shall assume that a ^ 0, 1, 2, — 
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3 Supersymmetrization of pairs of Hamiltonians 



Let us first introduce a compact notation and recollect the essence of the Witten's 
supersymmetric quantum mechanics. It introduces the auxiliary "conjugate" opera- 
tors A = d x +W and B = —d x + W in terms of the superpotential W. This factorizes 
certain Hamiltonians in full accordance with the old Schrodinger's recipe 0. For the 
two partner operators 



H (L) =p 2 + W 2 - W, H (R) =f + W 2 + W 



(5) 



it is easy to confirm that Hn\ = B ■ A while Him = A - B. In the Witten's language, 
the related two by two super-Hamiltonian and the two supercharges, 
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generate a representation of the superalgebra sl(l/l). The validity of the mixed 
anticommutation and commutation relations 



{Q,Q} = H, {Q,Q} = {Q,Q} = 0, 



[H, Q] = [H, Q] = 



mimics a supersymmetry between the bosonic and fermionic sectors of the Hilbert 
space. As already mentioned, the most appealing Fock-space picture of this type 
of the fermion - boson correspondence can be obtained in the "linear model" using 
the Hermitian one- dimensional harmonic oscillators. The scheme collapses after one 
moves to any dimension D > 1 {cf. the counterexample (488) in review 0}. 

Our £ / regularization of Schrodinger equations makes the spiked harmonic 
oscillators eligible for supersymmetric treatment. Let us return to our equation (01) 
and to its well known solvability in terms of Laguerre polynomials, 



C(r)=C<P( 



nl 



T(n + p + l) 



,13+1/2 



exp(-r 2 /2) -L<fV 



P = ±a 



(6) 
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A few algebraic consequences follow from the analytic VT regularity at e 7^ 0. We 
may pick up any real parameter 7 7^ 0, ±1, ±2, . . . and construct the superpotential 

.(7), 



-0 

The supersymmetric recipe gives the 7— numbered partner Hamiltonians (|5|). After 
the appropriate insertions they may be given the explicit and compact harmonic 
oscillator form, 

H™ = - 2 7 - 2, ifgj 27, a = | 7 |, 5=| 7 + 1|. (8) 

At all the non- integer 7 7^ 1/2 the direct computations reveal that the superscripted 
operators 

A M = d r + W M , B^ ] = -d r + W^\ 7 ^0,±1,... (9) 

act on our (normalized, spiked and VT— symmetrized) harmonic oscillator states in 
the transparent manner, 

A™ 4 7 i = Cl (n, 7 )4 7+1) , fh(n,i) = -2v^TT; 

B W 47+D = C2(n> 7) £ (7) i; C2(n> 7) = _2v^TT; 

A (7) £(-7) = Cs(n; 7) £ (-7-D ; C3(n; 7) = 2v /r7^; (10) 

B (7) 4" 7_1) = C4(n, 7) 4" 7) c 4 (n, 7) = 2^^- 

The former two rules were sufficient to define the well known one-dimensional anni- 
hilation and creation at a = 1/2. The latter two lines have to be added in order to 
move us to any a 7^ 1/2. Our operators mediate the supersymmetric mapping and 
import an explicit 7— dependence in C3 and C4. Its definitely most embarrasing conse- 
quence is the re-curring singularity at 7 =integer. This reflects the above-mentionend 
(unavoided) crossing of the energy levels. 
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4 Transitions between the different a 



We observe that each of the two partner Hamiltonians generates the two quasi- 
parity subsets of the equidistant energy values. At each n = 0, 1, . . ., the quadruplet 
of energies 

E$=4n, E$ = 4n-4ry, E$ = An + 4, E$=4n-4ry (11) 

corresponds to the respective wave functions 

Mi) r(-i) /•(7+1) /"(-7-i) 

The ordering of these states depends on the sign of the parameter 7 in a way sum- 
marized in Table 2. In this scheme, a mere shift of the energy spectrum (i.e., a 
coincidence of the Hamiltonians with a and a) is, obviously, exceptional. The re- 
lated parameter 7 e is given by the algebraic equation \ j e \ = \ j e + 1| with the unique 
solution 7 e = — 1/2. In the Hermitian limit e^Owe reproduce the one-dimensional 
example. Its supersymmetry is unbroken. 

A few nontrivial a ^ 1/2 examples of our generalized supersymmetric partner- 
ship are displayed in the pairs of the neighboring columns in Table 3. Its part 
(a) samples the two-way correspondence between the two different Hamiltonians 
H {1 ^ and H^ 3 / 2 \ The 7 = -3/2 VT supersymmetry between H {L) = H {3 ^ + 1 
and H( R ) = ift 1 / 2 ) + 3 is followed by the 7 = 1/2 correspondence between the 
doublet H (L) = - 3 and H {k) = -I. As a net result we obtain the ap- 

propriate generalization of the creation/annihilation pattern for the p— wave column 
of Table 1. 

At 7 = —1/2 we encounter the "degenerate" (and, in the present context, utterly 
exceptional) textbook 7 = —1/2 pattern 

a ■ dUSix) = V2n-^l£t?\x), a ■ C^ 2 \x) = -y/toi/M? (x) 

at CtT\x) = V2n-^TcS\x), at C^(x) = -V^C^ 2 \x) 
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where a ~ ^M/ 2 ) an d a" 1 " ~ B^ 1 / 2 ^. In a slightly re-ordered form, Table 3 (a) of- 
fers an alternative. Via the non- Hermit ian detour and limit e — > 0, another explicit 
annihilation pattern is obtained for the same s— wave oscillator. The new two-step 
mapping starts from the Hamiltonian = f/^ 1 / 2 ) — 3 and gives, firstly, its VT 
symmetrically regularized non-Hermitian supersymmetric partner = i/( 3 / 2 ) — 1 
[cf the last two columns in Table 3 (a)]. In the subsequent step, the similar partner- 
ship of the re-shifted H,^ = i/( 3 / 2 ) + 1 returns us to the original H,™ = ift 1 / 2 ) + 3 
[note that the shift differs from the one used in the first two columns in Table 3 (a)]. 



5 Creation and annihilation operators 

At any 7 7^ —1/2, our new supersymmetric pattern differs quite significantly from 
the Hermitian one. Firstly, in a fairly unusual way, the states can easily appear at 
a negative value of the energy [cf. the E^ = = —2 right-hand side sample 

in Table 3 (a) etc]. Secondly, in a way reflecting the same possibility, the usual 
"unmatching" (i.e., absence of the (^—subscripted partner) at the vanishing energy 
En?) = can occur for the excited states [cf., e.g., C^ 2 ^ — > in Table 3 (b); the 
existence of such an anomalous "excited vacuum" is particularly important for the 
non-equidistant spectra]. Thirdly, in a sharp contrast to the Hermitian models, the 
lowest state exists (i.e., remains normalizable) in both the spectra of and 
[cf., e.g., the 7 = 3/2 and E = -6 example C { Q ~ 3/2) -> C { ~ 5/2) in Table 3 (b)]. 

What can we expect from moving to the larger semi-integers a = Just 
a strengthening of the tendencies which were revealed in Table 3. Their features 
can be simply extrapolated. Thus, the large— 7 modifications of Table 3 will con- 
tain more lines at the bottom. For example, the 7 = 5/2 supersymmetry be- 
tween = ij( 5 / 2 ) — 7 and H( R ) = i/( 7 / 2 ) — 5 will induce a ground-state mapping 
£q~ 5 ^^—>£q~ 7 / 2 \ It appears at = — 10, witnessing just a continuing down- 

ward shift of the levels with the negative and decreasing superscripts. 
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Very similar conclusions can be drawn from the results generated at any real su- 
perscript value 7. We are near a climax of our study. A nice and transparent structure 
emerges form all the above PT supersymmetric assignments. Having spotted the 
difference between the regular and quasi-singular supersymmetries with a — 1/2 and 
a 7^ 1/2, respectively, we are ready to study more deeply the non-equidistant spectra 
corresponding to the non-integer values of 2a. 

Unless we reach the points of degeneracy a =integer, all our formulae remain 
applicable. The annihilation operators and their creation partners acquire the fac- 
torized, second-order differential form 

^(-7-1) . A (i) = 4(7-1) . A (-i) = A(a), 
B {-i) . 5(7-1) = B h) . s (-7-i) = A f (ai). 

At any a 7^ 0, 1, 2, . . ., they enable us to move along the spectrum of any harmonic 
oscillator Hamiltonian H^ a \ We get 

A(a). 4^i = c 5 (7i,/3)£f, 

At( a ).£(f) = c 5 (n,/3)4?i, 
c 5 (n, /3) = -4^/(n + l)(n + ^+l), /3 = ±a. 
This action is elementary and transparent. 

6 Summary 

We achieved a unified description of the spiked harmonic oscillators within the 
VT sy metric framework. 

• The general VT supersymmetric partnership has been shown mediated by the 
"shape-invariance" operators and 
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• At any non-integer a > the role of the general creation and annihilation 
operators for a given, single Hamiltonian has been shown played by their 
a— dependent and (3— preserving products A\a) and A(a), respectively. 

The only case where, up to a constant shift, the VT supersymmetric partners co- 
incide corresponds to the case where the poles in and vanish. This is the 
only case tractable without the use of the VT regularization. In this sense, the tra- 
ditional creation and annihilation using a^ = B(—l/2)/y2 and a = A{— l/2)/v / 2 is 
slightly misleading since these operators change the quasi-parity (3 to —(3. Such a 
correpsondence is not transferable to any non-equidistant spectrum with a ^ — 1/2. 

Our "natural" operators of creation A' (a) and annihilation A(a) are smooth near 
a = 1/2. Their marginal (though practically relevant) merit lies in their reducibil- 
ity to their regular (and hence, of course, state-dependent) first-order differential 
representation 

A(a) • 4f ) = (2r<9 r + 2 r 2 - An - 2(3 - 1) • 

A\a) ■ = (-2rd r + 2r 2 - 4n - 2(3 - 3) • Cf. 

Let us notice that the change of variables r — > y giving a simpler differentiation 
2rd r — > d y reads r = exp 2y and would result in the so called Morse form of the 
Hamiltonian. This could, in principle, indicate that the Morse potentials with VT 



symmetry [17 would deserve more attention. 
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Tables 



Table 1. 

A "disappearance" of wave functions in the Hermitian limit e — > at D = 3. 



harmonic oscillator 


energy 


linear (Z) = 1) 


spherical (D = 3) 

s — wave p — wave 




|4> = 4~ 1/2) 

|3) = C? /2) 

|2> = 4~ 1/2) 
li) = 4 1/2) 

|o) = 4~ 1/2) 


|1> = y/2C? /2) 
|0) = v/24 1/2) 


|1) = v^4 3/2) 

|0) = V2Cf /2) 


9 
7 
5 
3 
1 
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Table 2. 

Supersymmetry of harmonic oscillators at non-integer a — \^\. Superpotential (|7|) 
gives the partner Hamiltonians (|^) and quadruplets of energies E a < E\, < E c < E& 
at each n — 0, 1, 



range of a 


(0,oo) 


(0,1) 


(l,oc) 


parameters 
a = 7 + 1 

7 


a + 1 

a 


1 — a 
—a 


a — 1 
—a 


Hamiltonians 

rr(7) 
H (L) 

rrwJ 


iJ( Q ) - 2a 
#( fi ) - 2a 


- 25 
+ 2a 


+ 2a 
+ 2a 


energies 

= [of 4 a >] 

£ c = £fa [of 4")] 

e 6 = s (Ji) [of 4~«)] 
E a = [of 4-°)] 


4n + 4 

4n 
4n - 4a 
An - 4a 


An + A 
4n + 4a 
An + 4a 
4n 


4n + 4a 
An + 4a 
An + A 
An 
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Table 3. 

Supersymetry for the singular Hamiltonian if ( 3 / 2 ) = p 2 + (x — is) 2 + 2/(x — 
(a) 7 = -3/2 and -7 - 1 = 1/2. 



E, 



(L/R) 

14 

12 

10 

8 

6 

4 

2 





l"(L)) 



-(3/2) 
^2 

-(-3/2) 
^3 

^(3/2) 
£(-3/2) 



^(-3/2) 



^0 



(3/2) 



(-3/2) 



^0 



(-3/2) 



A(!/2) 



^2 

/•(-1/2) 
^2 

£(1/2) 
£(-1/2) 
r (l/2) 

r {-m 
^0 



£(3/2) 

-(-3/2) 
^2 

-(3/2) 
^0 

(-3/2) 



A-3/2) 
^0 



(b) 7 = 3/2 and -7 - 1 = -5/2. 



\n 



(L)l 



A( 3 / 2 ) 



n 



(«)/ 



n 



A<- 5 /2) 



(L)/ 



(«)/ 



8 
6 
4 
2 


-2 
-4 
-6 



-(3/2) 
^2 

A-3/2) 
*~3 

£(3/2) 

A-3/2) 
*~2 

-(3/2) 
^0 

£(-3/2) 



A-3/2) 



£(5/2) 

-(-5/2) 
^3 

-(5/2) 
^0 

-(-5/2) 
*~2 



£(-5/2) 



-(-5/2) 
^0 



^(3/2) 

-(-3/2) 
>-2 

-(3/2) 
^0 

£(-3/2) 



,(-3/2) 
-0 
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